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ABSTRACT 

Since four-fermion operators in strongly coupled QED are nonperturba- 
tively renormalizable, we analyze here the phase diagram and critical be- 
haviour of the Gauged Nambu-Jona Lasinio model. Our mean field approx- 
imation relates the critical exponents along the continuous phase transition 
line with the mass dependence of the chiral condensate in the Coulomb phase 
of standard noncompact QED. The numerical results for noncompact QED 
strongly suggest non mean field exponents along the critical line. 
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The possible existence of a non gaussian fixed point in strongly coupled 
non compact QED has been subject of extensive research during the last 
few years. It is well known that it has been a controversial subject for a long 
time. The present status of the field can be summarized in the following two 
points: 

1 . Triviality as in a logarithmically improved scalar mean field theory has 
been disproved by the numerical results 0. 

2. Unfortunately triviality a la Nambu-Jona Lasinio, which is the most 
natural way to analyze triviality in a theory with strongly coupled fermion 
fields, can not be disproved by the present day numerical data and it seems 
very difficult to get an important improvement in the quality of these data 
in near future. 

Therefore other alternative ways and new ideas are necessary in order to 
get some progress in this field. This is the reason why the gauged Nambu- 
Jona Lasinio (GNJL) model has become increasingly interesting in recent 
time. In fact if a non gaussian fixed point exists in non compact QED, the 
naive dimensional analysis does not applies. Therefore operators of dimen- 
sion higher than four, which are non renormalizable in perturbation theory, 
could acquire anomalous dimensions and become renormalizable [El. 

The lattice action for the GNJL model with noncompact gauge fields 
and staggered fermions reads 



S= 2 £ ®/U n ) + X A (0)X + m XX ~ G^XnXnXn+pXn+v (1) 
n,fi<u n,fi 

In the chiral limit, m = 0, this action is invariant under the continuous 
transformations 

Xn - Xn e*-» ni+ - +nd Xn - Xne^-^ + - +nd (2) 

which define a continuous chiral U(l) symmetry group. 

The main technical difficulty when computing vacuum averages in the 
GNJL model comes from the fact that the action (1) is not a bilinear of 
the fermion fields. The standard procedure consists in the introduction of 
an auxiliary vector field which allows to bilinearize the fermion action. The 
prize to pay for that is that we have one more field to include in the nu- 
merical simulations of this model that besides the number of free parameters 
(j3, m, G), makes it difficult to analyse this model with reasonable computer 
resources [II. 
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As a first approach, we can do a standard mean field approximation which 
has the advantage of bilinearizing the action (1). Following the mean field 
technique, we do in (1) the following substitution 

G^XnXnXn +t xXn+n ~> 2rfG( XX ) XnXn (3) 
n,fi n 

where d is the space-time dimension. The action (1) becomes in this way 
a bilinear in the fermion fields and the path integral over the Grassmann 
variables can be done by mean of the Mathews-Salam formula. 

The v.e.v. of the chiral condensate after the substitution of the mean field 
approximation (3) in the action (1) is given by 



( _ /He-lEe^(n) de t[ A + ( m -8G(xx))/]^r A+(m ^ G( _ x))J 

J[d9}e-^ e2 ^det[A + (m-8G(xx))I] 



In the chiral limit m — 0, equation (4) becomes 

■ 1 v/2 1 



(4) 



(XX) = 16G(xx) ( ttE , 2 , RAn2l - ,2 ( 5 ) 
\ V j=1 Xj + 64G 2 (xx) I 

where the sum in (5) runs over all positive eigenvalues of the massless Dirac 
operator and the integration measure in the v.e.v. includes the fermionic 
determinant of standard noncompact QED, evaluated at the effective mass 
m = -8G(xx)- 



1. The phase diagram 

Equation (5) is always verified if ( XX ) = 0> an d this is the only solution 
in the symmetric phase. In the broken phase where (xx) 7^ 0> the v.e.v. of 
the chiral condensate will be given by the solution of the following equation 

/ 1 v/2 1 \ 
l = 16G(-y ? ) (6) 

which gives for the critical line, where the chiral condensate vanishes contin- 
uously, the following expression 

GJP) = ^ (7) 

V \ ^3=1 ~>? / 
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The existence of this critical line in the GNJL model was discovered 
some time ago Q in the continuum formulation using the quenched-ladder 
approximation. 

In Fig. 1 we present our numerical results for the phase diagram in the 
(3, G plane. The critical line has been obtained by computing numerically 
the v.e.v. of the sum of the inverse square eigenvalues (eq. (7)), which is 
proportional to the chiral transverse susceptibility of the standard noncom- 
pact QED in the chiral limit. The numerical simulations where done using 
the MFA approach ||, which allows to do computations in the chiral limit. 
We refer the interested reader to the extended bibliography on this subject 
[|] and specially to the ref. where the computation of the chiral suscepti- 
bility and the determination of the critical coupling in noncompact QED is 
discussed in detail. 

2. The critical exponents 

The phase diagram of Fig. 1 is in good qualitative agreement with the 
corresponding phase diagram obtained in the quenched-ladder approximation 
[|]]. Using this analytical approach, a line of critical points with continuously 
varying critical exponents was found in |7J], the intersection point of this line 
with the G = axis corresponding to an essential singularity 0. 

Later on, numerical simulations of noncompact QED disproved the essen- 
tial singularity behavior, putting in evidence the limitations of the quenched- 
ladder approximation. Since our approach contains less approximations, we 
do hope to get more reliable results for the critical exponents. 

In order to extract the critical exponents, we will start from the key equa- 
tion of state relating the order parameter with the external magnetic field m 
and the gauge and four fermion couplings. Using the previous notation we 
can write 

{xx) = -2mF{P,m) (8) 

where the right hand side in (8) is just the chiral condensate in full non- 
compact QED evaluated at the gauge coupling value (3 and fermion mass 
m. Concerning critical exponents the interesting physical region, as follows 
from the phase diagram of Fig. 1, is fi > (3® (Coulomb phase of noncompact 
QED). 

Since we are interested in the critical region (m — > 0, (xx) ~ * 0), we will 
analyze the behavior of F((3,fh) in the fh — > limit. In this limit we can 
write 
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F(0,m) = F(0,O) + Bm UJ + 



(9) 



The second term in (9) can contain also logarithmic contributions and 
F((3, 0) is one half of the massless transverse susceptibility in noncompact 
QED. Equation (9), after the substitution of m by m - 8G(xx), implies 
the following behavior for the chiral condensate in the m — > limit 

(xx)^m^TT (io) 
and therefore the u and 5 exponents are related by the equation 

5 = u + l (11) 

A straightforward calculation allow to compute also the magnetic f3 m and 
susceptibility 7 exponents, the final result being 

Pm = ~, 7=1 (12) 

The hyperscaling relation 7 = j3 m (6 — 1) is verified, as follows from (12). 

The determination of the order parameter critical exponents in our mean 
field approach reduce therefore to the determination of the u exponent which 
controls the mass dependence of the chiral condensate in the Coulomb phase 
of noncompact QED. In the (3 — ► 00 limit of noncompact QED, the theory is 
free and the chiral condensate can be analytically computed. The well known 
result in this case (u — 2 plus logarithmic corrections) implies mean field 
exponents for the end point of the phase transition line, with the following 
behavior for ( XX )/3=oo,g=g- 

m ~ (xxflog(xx) (13) 

In the general case, the chiral condensate in the Coulomb phase of non- 
compact QED ((xx)ncqed) can be parameterized as follows 

( XX )ncqed = A([3)m + B(f3)m" +1 + ... (14) 

The first contribution in (14) is linear in m, as follows from the fact 
that the massless transverse susceptibility is finite in the Coulomb phase of 
noncompact QED. The next contribution can have logarithmic corrections, 
as happens in the (3 — > 00 limit where it becomes m 3 log m. In order to extract 
the to exponent from the numerical simulations, we can use the results for 
the massless chiral transverse susceptibility to fix A(/3) in (14) and fit 
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the numerical results with eq. (14). This procedure has the inconvenient 
that higher order contributions in (14) can induce systematic errors in the 
determination of uj. 

A better strategy exploiting the potentialities of the MFA method is 
to compute vacuum expectation values of operators which can be consid- 
ered as generalizations of one of the contributions to the massless nonlinear 
susceptibility. More precisely we have defined \ q by the expression 




When q = 4, we get one of the contributions to the standard massless 
nonlinear susceptibility. In the general case we can write this vacuum ex- 
pectation value as an integral over the spectral density of eigenvalues in the 
following way 




and if the density of eigenvalues p(A) behaves like \ p near the origin, \q will 
diverge when q > p + 1. In such a case and for lattices of finite size, we 
expect for \q the following type of divergency with the lattice size L 

Xq „ L a( q -p-l) (17) 

where a in (17) is some positive number. 

An interesting thing to notice is that the p-exponent which controls the 
small A behavior of the spectral density p(X), can be related with the uj 
exponent by the following equations 

co = p — l{p < 3) 

uj = 2(p > 3) (18) 

These relations allow to extract the uj exponent from the finite size be- 
havior of the generalized nonlinear susceptibility x q . 

In Fig. 2 we have plotted our results for the inverse of the generalized 
nonlinear susceptibility \q against the inverse lattice size for q- values running 
from 2 to 4 and (3 = 0.237. The solid lines in this figure correspond to a fit 
of all the points at any fixed q with the function 
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X -\L) = a q + b q L- c (19) 

The results reported in this figure show how the infinite volume limit 
of the inverse generalized nonlinear susceptibility vanishes at large q and is 
different from zero at small q, as expected. Fig. 3 is a plot of the extrapolated 
values of Xq 1 (thermodynamical limit) against q. The critical value of q 
at which x~ x vanishes can be estimated from these results. Hence we get 
q c ~ 2.5 at j3 — 0.237, which implies p ~ 1.5 and to ~ 0.5. Using now 
the relations (11), (12) the following results for the order parameter critical 
exponents can be derived 

5-1.5, f3 m ~2, 7=1, (20) 

values which are clearly outside the range of the mean field exponents. 

Even if at first sight non mean field exponents in a mean field approach 
seems to be a paradox, there is no contradiction between both statements. 
In fact we have applied the mean field approximation to the fermion field 
but fluctuations of the gauge field are taken into account in our numerical 
simulations. In the infinite (3 limit, where the gauge field is frozen to the free 
field configuration, we get mean field exponents. However fluctuations of the 
gauge field at finite (3 seem to play a fundamental role in driving critical 
exponents to non mean field values. 

The picture which emerges from this calculation is that the critical expo- 
nents change continuously along the critical line of Fig. 1 from their mean 
field values (end point of the critical line) to some non mean field values at 
the critical point of noncompact QED. The 5 exponent approaches its mean 
field value (5 = 3) from below whereas the magnetic exponent approaches 
its mean field value (/3 m = 0.5) from above |§. Our results for several values 
of the gauge coupling f3 suggest also that the value of 5 increases systemati- 
cally along the critical line with increasing /3, in contrast with the magnetic 
exponent results which are systematically decreasing with (3. 

In conclusion we do believe our qualitative picture is realistic. In fact it 
is hard to think that non mean field exponents in a mean field approach will 
become mean field exponents after removing the mean field approach. 
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Figure captions 



Figure 1. Phase diagram of the GNJL model in the (3, G plane. 

Figure 2. Logarithm of the nonlinear susceptibility against the logarithm 
of the lattice size for lattice sizes 4,6,8, 10 and j3 = 0.237. 

Figure 3. Inverse generalized nonlinear susceptibility against the inverse 
lattice volume at (3 — 0.237. 

Figure 4. Infinite volume limit of the generalized nonlinear susceptibility 
against q at ft — 0.237. 
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